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We analyze the dynamics of model supercooled liquids in a temperature regime where 
predictions of mode coupling theory (MCT) are known to be valid qualitatively. In 
this regime, the Adam-Gibbs (AG) relation, based on an activation picture of dynam¬ 
ics also describes the dynamics satisfactorily, and we explore the mutual consistency 
and interrelation of these descriptions. Although entropy and dynamics are related 
via phenomenological theories, the connection between MGT and entropy has not 
been argued for. In this work we explore this connection and provide a microscopic 
derivation of the phenomenological Rosenfeld theory. At low temperatures the over¬ 
lap between MGT power law regime and AG relation implies that the AG relation 
predicts an avoided divergence at Tc, the origin of which is traced back to the van¬ 
ishing of pair conhgurational entropy, which we hnd occurs at the same temperature. 
We also show that the residual multiparticle entropy plays an important role in de¬ 
scribing the relaxation time. 
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I. INTRODUCTION 


In the study of liquid state physics, the structure of the liquid, which is often described 
primarily by the two body radial distribution function (rdf), has always played a central 
role. The structure can not only describe the thermodynamic properties of the liquid like 
the energy and pressure, under certain theoretical frameworks like the mode coupling theory 
the structure can also determine the dynamicsii^. In a series of paper Berthier and Tarjus 
have described the behaviour of two systems with different interaction potentials, namely, 
the Lennard-Jones (LJ) and the Weeks-Chandler-Andersen (WCA) potentials. Although 
at the same temperature and density the structures of these systems are very close, the 
dynamics display signihcant differences at low temperatures^^-. These studies questioned 
the role of structure in determining the dynamics. Coslovich has shown that although the two 
body radial distribution function of these two systems are quite similar, triplet correlations 
are signihcantly different^. He has also shown that the LJ system has more pronounced 
local ordering^. In supercooled liquids these locally preferred structures are known to form 
correlated domains which have been argued to give rise to the slow dynamics^. An estimation 
of this length scale of the domains and its connection to the relaxation timescale is a topic 
of ongoing research^^J^. One such study by Hocky et al. has shown that the point-to-set 
correlation length of the LJ system is larger compared to that of the WCA system and that 
this difference in correlation length can account for the difference in dynamics of the two 
systemsl^. From these studies one may conclude that the difference in dynamics primarily 
comes from many body correlations. However, in a recent study by some of us it has been 
shown that two body correlation information is good enough to capture the difference in the 
dynamics between the two systems. The study also reveals that the divergence temperature 
at which an approximation to the conhgurational entropy using pair correlation alone goes 
to zero, is similar to the mode coupling theory (MCT) transition temperature, T^. As 
mentioned before, MCT is a microscopic theory where the structural inputs determine the 
dynamics. Although entropy and dynamics are related via phenomenological Rosenfeld^ 
and Adam-Gibbs (AG)^ relations at high and low temperatures respectively, MCT does not 
have any apparent connection to entropy. Thus it is of great interest to try to understand 
the origin of the coincidence of the MCT divergence temperature and the temperature where 
pair conhgurational entropy goes to zero. 


2 


At normal liquid temperatures, a semi quantitative correlation between the dynamics 
(transport properties) and thermodynamics (excess entropy), proposed by Rosenfeldl^'i^, 
has been extensively studied in recent times, where the relaxation time r can be written as, 

r(T) = Cexp [-KSex] (1) 

Here C and K are the constants. Since the pair entropy S 2 , which is obtained only from the 
pair correlation function, accounts for 80% —90% of the excess entropy^i^^, many simulation 
studies have replaced S^x by S 2 and have shown that even with S 2 the transport coefficients 
follow Rosenfeld scaling^'^^^— . 

Bagchi and coworkers used Zwanzig’s rugged energy landscape model of diffusion^^ and 
by connecting the ruggedness to the excess entropy have provided a derivation of Rosen¬ 
feld relatioii^^. Samanta et am have shown that under certain approximations the diffusion 
coefficient as obtained from MCT follows Rosenfeld scaling. Das and coworkers have per¬ 
formed microscopic MCT calculations which show that the diffusion values thus obtained 
can be htted to Rosenfeld scaling^^. Some of these studies have reported that the scaling 
parameter is not unique, hence the whole temperature region cannot be htted to a single 
straight line^^i^^. 

Although Rosenfeld scaling holds at high temperature, it is known to breakdown at low 
temperatures even with multiple scaling parameters. At low temperatures the correlation 
between the transport coefficients and entropy is usually described by the well known Adam- 
Gibbs relation!^, 

r(r) = r^exp , (2) 

where Sc is the conhgurational entropy of the system. For a wide range of systems, the AG 
relation is found to hold^‘^‘2^ below a moderately high temperature referred to as the onset 
temperature of slow dynamics. 

In this paper, we explore the connection between dynamics, characterization of structure 
as contained in the pair and higher order correlations, and entropy, and relations between 
them as described by MCT and the AG relation, using computer simulations of two model 
liquids and analytical results that seek to relate descriptions of dynamics in terms of struc¬ 
ture, and entropy. Our present study shows that the AG theory, which is based on activation 
dynamics can completely describe the mode coupling theory (MCT) power law behavior in 
the region where the latter is found to be valid. An earlier study also observing similar 
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overlap regioi>2£ argued that the observation supports the hypothesis that a direct relation 
exists between the number of basins and their connect ivity22i2i. In this work to understand 
the above mentioned observations, we explore the connection between mode coupling theory 
(MCT) and entropy and discuss different predictions of MCT in the light of entropy. We 
also analyze the different roles of pair and many body correlations. 

Although MCT makes predictions about dynamics in both Rosenfeld and AG temperature 
regimes, no connection between MCT and entropy has been argued for, except for one 
study 2 ^, as mentioned earlier. We show that, under some assumptions, the memory function 
in the MCT equation for structural relaxation is related to the pair excess entropy, thus 
providing a microscopic derivation of the phenomenological Rosenfeld expression for the 
structural relaxation time, r. Our study also can explain the origin of the temperature 
dependence of the Rosenfeld parameter. The origin of higher relaxation time and higher 
activation energy as predicted by MCT is also obtained from the analysis of the memory 
function. 

As mentioned above the AG expression for relaxation times and the MCT power law form 
overlap in a certain temperature regime. The AG relation is valid for a wide temperature 
range which includes the range in which the MCT power law prediction holds. Thus in the 
MCT regime, the relaxation time follows both the AG and MCT behaviour. Our study 
reveals that the origin of the avoided divergence like behaviour (as given by MCT power 
law) in the AG relation is related to the vanishing of the pair conhgurational entropy. 
However we show that the pair conhgurational entropy, although predicting the correct 
MCT transition temperature, by itself cannot predict the MCT power law behaviour. The 
residual multiparticle entropy (RMPE) plays an important role in providing the correct 
temperature dependence of relaxation times. We also hnd a connection between the AG 
coefhcient (A), pair thermodynamic fragility {Kt2) and MCT critical exponent (7). We 
show that although both ‘A’ and Kt2 are dependent on density, their ratio which is related 
to 7 is density-independent. 

The paper is organized as follows: The simulation details are given in Sec. II. In Sec. Ill 
we describe the methods used for evaluating the various quantities of interest and provide 
other necessary background. In Sec-IV we report some observations that motivate our 
analytical results which are described in Sec-V. In Sec-VI we present additional numerical 
results and their analysis. Sec. VII contains a discussion of presented results and conclusions. 
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II. SIMULATION DETAILS 


We have performed molecular dynamics simulations of the Kob-Andersen model which 
is a binary mixture (80:20) of Lennard-Jones (LJ) particles and the corresponding WCA 
versioii^*^. The interatomic pair potential between species i and j, with i,j = A,B, Uij{r) 
is described by a shifted and truncated Lennard-Jones (LJ) potential, as given by: 


U„ir) = 


0 , 


(c) \ , (c) 

r)Aaij,eij), r<r\/ 




n 


(3) 


where aij, eij) = 4ejj[(crjj/r)^^ —(cTjj/r)®] and = 2.5(Tjj for the LJ systems and 

is equal to the position of the minimum of for the WCA systems. Length, temperature 
and time are given in units of an, ksT/en and r = yjimiah/en), respectively. Here we 
have simulated Kob Andersen Model with the interaction parameters an = 1.0, ai 2 =0.8 
,a22 =0.88, en =1, =1.5, €22 =0.5, mi = m2=1.0 . 

The molecular dynamics (MD) simulations have been carried out using the LAMMPS 
package^. We have performed MD simulations in the canonical ensemble (NVT) using NosA 
Hoover thermostat with integration timestep 0.005r. The time constants for NosAHoover 
thermostat are taken to be 100 timesteps. The sample is kept in a cubic box with periodic 
boundary condition. System size is A = 500, Na = 400 (N = total number of particles, Na 
= number of particles of type A) and we have studied a broad range of density p from 1.2 
to 1.6. For all state points, three to hve independent samples with run lengths > IOOtq, (r^ 
is the a- relaxation time) are analyzed. 
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III. DEFINITIONS AND BACKGROUND 

A. Relaxation time 


We have calculated the relaxation times from the decay of the overlap function q(t), using 
q{t = Ta,T)/N = 1/e. The overlap function is dehned as 


m) ^ 

N N 


drp{r, to)p{r,t + to) 


- ri{t + to)) 

i=l j=l 


N 




. 2=1 


( 4 ) 


i j^i 

The overlap function is a two-point time correlation function of local density p{r,t). It has 
been used in many recent studies of slow relaxatiort^l. In this work, we consider only the 
self-part of the total overlap function (i.e. neglect the i ^ j terms in the double summation). 
This approximation has been shown to be a good approximation to the full overlap function. 
So, the self part of the overlap function can be written as. 


N 


{(lit)) ~ ( ^d{ri{to) - ri{t + to)) 


( 5 ) 


, 1=1 


The 6 function is approximated by a window function u}{x) which dehnes the condition 
of overlap between two particle positions separated by a time interval t: 


N 


{(lit)) ~ I) 

\ i=l I 

uj{x) = l,x < & implying overlap 
= 0, otherwise 


( 6 ) 


The time dependent overlap function thus depends on the choice of the cut-off parameter 
a, which we choose to be 0.3. This parameter is chosen such that particle positions separated 
due to small amplitude vibrational motion are treated as the same, or that is comparable 
to the value of the MSD in the plateau between the ballistic and diffusive regimes. 

Relaxation times obtained from the decay of the self intermediate scattering function 
Fs{k,t) using the dehnition Fs{k,t = Ta,T) = 1/e at /c ~ 27i/rmax, where r^ax is the hrst 
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maximum of the radial distribution function . The self intermediate scattering function is 
calculated from the simulated trajectory as 

N 


Fs{k,t) = — ( ^exp(-ik.(ri(f) -ri(0))) 


(7) 


. 2=1 


Since relaxation times from q{t) and Fs{k,t) behave very similarly at low temperature 
we have used the time scale obtained from q{t). As q{t) cannot be calculated from MCT, at 
high temperatures where we compare the simulation result with the relaxation time obtained 
analytically from MCT we have computed relaxation time using Fs{k,t). 

B. Static Structure Factor 


We measure the partial structure factor Sa^{k) which are needed as input for the MCT 
calculations. They are dehned as 


Saf^ik) = 


\/NaN/3 


5^5^exp(-*k.(r“-rf)) 


( 8 ) 


C. Mode coupling Theory 

Many properties of a glass forming liquids can be explained by the well known mode 
coupling theory of the glass transition (MCT). This microscopic theory can give a qualitative 
description of dynamical properties (such as temperature dependence of relaxation time) if 
the static structure of the liquid is known and many experiments and simulation results has 
shown that MCT predictions hold good in the temperature regime of initial slow down of 
dynamics^. The equation for the intermediate scattering function (j){k, t) is given by 


0(/c, f) + r0(/c, f) + fi^0(/c, f) + J dt'Ai{t — t')^{k,t') = 0 

where and memory function of (j){k,t) can be written as : 


(9) 


M{k,t) = 


dq 


■V'fc (q,k-q)^(fc)^(| k-q \)S{q)(p{q,t)(p{\ k-q |,f) (10) 


2pk‘^ J (27r)^ 

where k — q = p and 14(q, p) = [k.qpC'(g) + k.ppC(p)] . 

For the self intermediate scattering function a similar equation may be written, as 

0s(/c, t) + r0s(/c, t) + Ql(j)s{k, f) + fig [ dt'Msit - t')(f)s{k, t') = 0 


( 11 ) 
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where fig = and memory function Ais{k,t) can be written as 

We need the static structure factor to solve these equations, which is obtained from 
computer simulation, by EqJHl The temperature dependence of the system enters in MCT 
through S(q) and since we need very precise S(q) near the MCT transition, we simulated 
S(q) at three temperatures around the transition point and used them to create the structure 
factors at intermediate temperatures by quadratic interpolation method as described 

To solve the Eq|TT] we need (p{k,t) as an input, which can be taken by solving EqlHl 
Using these expressions we have calculated the relaxation time, tmct from the relaxation 
of (j){k, t) at 1 /e. 


D. Configurational Entropy 

Conhgurational entropy. Sc per particle, the measure of the number of distinct local 
energy minima, is calculated^ by subtracting from the total entropy of the system the 
vibrational component: Sc(T) = StotaiiT) — Smb(T)^’^. The total entropy of the liquid is 
obtained via thermodynamic integration from the ideal gas limit. Vibrational entropy is 
calculated by making a harmonic approximation to the potential energy about a given local 
minimum. 


E. Pair Configurational Entropy 


To get an estimate of the conhgurational entropy as predicted by the pair correlation we 
rewrite Sc in terms of the pair contribution to conhgurational entropy Sc 2 ^, 

Sc = Sid + Sex ~ Syib = Sid + 5*2 + AS* — S^ib = Sc2 + AS* (13) 


Where Sc 2 = Sid + S 2 — S^ib- Sex can be expanded in an inhnite series. Sex = -S '2 + S '3 +.= 

S 2 + AS* using Kirkwood’s factorizatiort^ of the N-particle distribution functiort^Ai. Sn 
is the “n” body contribution to the entropy. Thus the pair excess entropy is S 2 and the 
higher order contributions to excess entropy is given by the residual multiparticle entropy 
(RMPE), AS = Sex - S 2 





IV. OBSERVATIONS 


As the liquid is supercooled, the Rosenfeld scaling, observed to be valid at normal tem¬ 
peratures, is known to break down^^. However, in this regime the Adam-Gibbs relation is 
found to holdi^i^l. The Adam Gibbs relation explains the behaviour of dynamical property 
like relaxation time using conhgurational entropy which is a thermodynamical property. So 
this relation connects thermodynamics and dynamics for low temperature liquids. In the 
Adam-Gibbs relation it is not the excess entropy but the conhgurational entropy which 
dictates the dynamics. 



FIG. 1. The power law behaviour of relaxation time of numerieal simulation, t, as predieted by 
MCT (Eq\Tf\) appears as a straight line for a certain region — 1) < for both the 

systems at all densities. The critical exponent 7 is obtained from the slope of the linear fit. For 
elarity, data at different densities are vertieally shifted. 

Although microscopic MGT shows a divergence of the relaxation time, r, at a much 
higher temperature^ than the glass transition temperature, the power law behaviour of r 
as predicted by MGT is found to be valid in a range of low temperatures. Similar to the 
earlier studies^*^, the power law behaviour of simulated r we compute is well described by 
an algebraic divergence given by, 

T ~ (T - X)--’ ~ (T - i)--, ( 14 ) 

c 

For all the densities we study, as shown in Fig HI in a certain region of temperature, 
(10“^ < — 1) < 10°), the relaxation time ,r, for both LJ and WGA systems follow the 

MGT power law behaviour. On the other hand, the Adam Gibbs relation is also valid for 
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FIG. 2. plotted against — 1) and as predicted by Eall5\ the plot is linear in the region 

10“^ < (^ — 1) < 10® validating our claim that MCT power law region overlaps with AG region. 

all the systems in this region (0 < — 1) < 10®)^. Thus we hnd that the temperature 

range where MCT like behaviour is predicted completely overlaps with the range where 
Adam-Gibbs relation is found to be valid. As mentioned in the Introduction this overlap 
regime has earlier been reported for other systems^E. As in this temperature regime Tq can 
be described both by MCT power law behaviour and and by the AG relation we can write, 

^oc-7 ln(^-l) (15) 

In FigJ2]we show that is linear when plotted against ln(^ — 1) in the region 10“^ < 
(^ — 1) < 10® validating the statement that MCT like divergence region overlaps with AG 
region. 

Since the conhgurational entropy has a hnite value at the MCT transition temperature, 
Tc, the AG relation is not expected to predict a divergent relaxation time at this temperature. 
In order to investigate the origin of this avoided transition, we consider the separation of 
the conhgurational entropy into pair and many body parts as described earlier (sec-3. 5)^. 
We hnd that the temperature dependence of (*S'c 2 ) is given by (Figj3]), 


TS,2 = Kt2{^ - 1) (16) 

J-K2 

where Kt 2 is the pair thermodynamic fragility and Sc 2 vanishes at the Kauzmann temper¬ 
ature Tk 2 is obtained from the linear ht of TSc 2 vs T plot at Sc 2 = 0. As reported 
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FIG. 3. The temperature dependence of pair configurational entropy (Sc 2 ) to determine Kauzmann 
temperature Tk 2 - Tk 2 values are given in Table\^ 


earlier we find that for all the systems studied in this work the Kauzmann temperature for 


Sc 2 is very close in value to the MCT transition temperature (Tabled]). 


TABLE I. Tffi- and Tk 2 values are tabulated below. For all the systems studied here, the Kauzmann 
temperature for Sc 2 is quite similar to the MCT transition temperature. 



p = 1.2 

p= 1.4 

p = 1.6 


Tc 

Tk2 

Tc 

Tk2 

T 

Tk2 

LJ 

0.435 

0.445 

0.93 

0.929 

1.76 

1.757 

WCA 

0.28 

0.268 

0.81 

0.788 

1.69 

1.696 


Thus, although Sc is hnite at the estimated MCT T^, Sc 2 vanishes at Tk 2 which coincides 
with Tc- 

V. ANALYTICAL RESULTS 

Our study shows that the AG theory which is based on activation dynamics can com¬ 
pletely describe the mode coupling theory (MCT) power law behavior in the region where 
the latter is found to be valid (FigJ2|). However, the microscopic picture for Mode Coupling 
Theory (MCT) and the Adam Gibbs (AG) relation are different. Either from the heuristic 
arguments of Adam and Gibbs, or from the Random First Order Transition (RFOT) deriva¬ 
tion, the AG relation is obtained from an activation picture of the dynamics, whereas the 
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MCT does not correspond to activated dynamics. This leads to the question of the role of 
entropy in MCT which will be the focus of this section. 

A. Entropy and MCT 

In the fc —?■ 0 limit, the memory function, Ai{k,t) in EqJTOl can be rewritten as. 


= dqx q‘^{S{q)-if X (17) 

In the Schematic MCT the 0(g, t) is usually decoupled from q, as in the memory function, 
Ai{k,t), the dominant contribution comes from the hrst peak of S(q)^‘^. Here we consider 
similar decoupling, however do not restrict ourself to hrst peak of S(q). Thus we write EqITTi 
as 


M{k,t) = 


S{k) 


dq{S{q) - 1)^ 


4p/c(27r)3 [7o 

By writing S(q) in terms of g(r) we can rewrite EqJTSl as 


X (j)‘^{k, t) 


M{k,t) = X 27ip 


drr‘^[g{r) — 1)" 


(f{k,t) 


(18) 


(19) 


Replacing Ai{k,t) from Eq{T9] in Eql9] and considering over damped limit by omiting the 
explicit ‘k’ dependence of 4>{t), Eql9]can be written in schematic form as 

0 (t) + f2^0(t) + f dt'(f{t')^{t — T) = 0 (20) 

Jo 

Where we can identify the coupling parameter A from Eq{T9] as 


A = ^ X / drr\g{r) - 1)= = (21) 

2k J 2k k^ 

Where we call S 2 approx as the approximate pair entropy. The choice of calling it entropy will 
become clear in the next analysis. 

We note that the two body pair entropy is given by^^. 


ks 


drr‘^{g{r) ln 5 f(r) — [g{r) — 1]} 


Expanding the logarithmic term for g{r) > 0 we get 

q roo 1 

r^ = -2^p/ dTTfg{T)-lf 

ks Jo {g{r) + 1) 


+ H 


( 22 ) 


(23) 
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where in ‘H’ we put the higher order contributions. The Figj4] shows that the primary 
contribution comes from the first term of Eql23l 

In the above equation we note that r^[g{r) — 1]^ varies strongly compared to l/{g{r) + 1). 
In the later if we consider g{r) ^ 1 we can write 

Q poo Q 

= -2np / drr^[g{r) - 1]^ 2^ - 2H (24) 

Kb Jo Kb 

Our numerical analysis shows that for all the systems studied here S 2 approx vs S 2 is indeed 
linear (FigJS]) with a slope ~ 2.5. Thus the coupling constant A is related to the pair entropy, 



FIG. 4. The function (^ — H) of Eq\^ and 82 /kB are plotted as a function of temperature. The 
plot shows that the primary contribution comes from the first term of the expansion. 

^ X _ H) (25) 

where rris is the slope obtained from S 2 approx vs S 2 plot. 

The MCT relaxation time from schematic model^ is given by 

r ~ (1 - A)"^ (26) 

Note that the power law behaviour of relaxation time r (as given by Eqj26]) changes to expo¬ 
nential dependence of r under generalized MCT formalism^, when the coupling parameter 
is considered to be the same for all higher order terms and frequency G ~ 1. With these 
conditions r can be written as 

^ = ^(exp(A) - 1) ~ ~ C"exp(iF'^ 2 ) (27) 
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FIG. 5. S 2 approx is plotted against S 2 and it shows a linear behaviour with a slope « 2.5. 

The second equality is written by replacing A from Eqj25l Where C and K' are not a 
constants, rather have a temperature dependence. 

Earlier study of diffusioii^^ and our present microscopic derivation of the Rosenfeld rela¬ 
tion for relaxation time r shows that similar to Rosenfeld prediction, the MCT also predicts 
it to be an universal scaling law for all transport coefficients. 

VI. NUMERICAL RESULTS 
A. Rosenfeld scaling and MCT 

In this section we analyze the MCT results in the light of Rosenfeld relation. We hnd that 
the relaxation time as obtained from microscopic MCT, tmct when plotted against A does 
not follow the power law ((1 — A)”'’') or exp{X) dependence in the whole temperature region. 
Usually it is found^'^ that both tmct and r (relaxation time obtained from simulation) 
when plotted against S 2 does not show a single straight line. In Figj6] we plot the tmct 
calculated from Eq.9-12 against S 2 which shows two linear regimes. The origin of this break 
or the temperature dependence of the Rosenfeld parameter ‘JC’ is not known. 

Our analysis of EqJ25] shows that the Rosenfeld parameters are related to the static 
structure factor S{k). Thus the temperature dependence of S{k) leads to the temperature 
dependence of Rosenfeld parameter ‘JC’. However since S{k) changes continuously with 
temperature, it should lead to a similar temperature dependence of K'. That a continuously 
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FIG. 6. (a) The relaxation time obtained from microscopic MCT, tmct is plotted against 82 - The 
dashed lines illustrate the two different Rosenfeld regimes, (b) Plot of X vs. S 2 . This also shows 
two different linear regimes. For clarity r and A are shifted by 1.2, 2.8 and -0.2, -O .48 for the 
systems of p = 1.4 and p = 1.6 respectively. The break in the slope for both the plots are illustrated 
by vertical dash-dot lines. We show that for a fixed density the break for both tmct and A are at 
the same S 2 value. 


changing ‘ JG’ is not needed to describe the observed behaviour but two distinct values suffice 
can be seen when we plot A against S2 (Figj6]-b), where we see that there is a break in the 
slope and it happens at the same S2 value where r against S2 shows a break in slope. 

Next we show that the value of S 2 approx and its temperature dependence as compared 
to Sex can explain i) the larger values of tmct as compared to r— ii) the higher values of 
activation energy as predicted by MCT-. When and are obtained by fitting r 

and Tmct to Arrhenius expression fEq j28l) we find values shown in Table 2, and in Fig.7. 

r ~ To exp ^ (28) 

Fig|7] shows that at all densities for both the systems S2approx is smaller than Sex and has 
a much stronger temperature dependence. Using Rosenfeld Expression we can write 

T{T) = CeM-KSex) (29) 


Now if we replace Sex by S 2 approxi keeping C and K same, we get 


^2approx 


C GXp( KS2approx) 


(30) 


The C and K are obtained from linear fits of logarithmic of simulated relaxation time against 
excess entropy. Since S 2 approx « Sex, the study shows that T 2 approx » Similar to that 
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FIG. 7. Sex and S 2 approx are plotted as a function of temperature for LJ and WCA systems at 
densities 1.2, I .4 and 1.6. For all the systems S 2 approx has stronger temperature dependence and 
has smaller value than Sex- In inset we plot T 2 apporx as obtained from EalWi. It shows that T 2 approx 
has higher value and a larger slope leading to higher activation energy as compared to r. Activation 
energies are tabulated in TableW\ 

predicted by microscopic MCT (Eqlll [TT]), the Eq values for T 2 approx are higher, which are 
given in Table UTl 

Although the results obtained from T 2 approx shows the correct trend, it can not match the 
parameters as obtained from tmct- We note that the T 2 approx is a prediction obtained from 


16 


















































TABLE II. Eq are tabulated for different systems. We show that the Eq values are higher for MCT 
as well as for approximate calculation. As we can not calculate overlap function from MCT, for 
comparison of Eq values we use simulated Fs{k,t) 



p = 1.2 

p= 1.4 

p = 1.6 


LJ 

WCA 

LJ 

WCA 

LJ WCA 

Tpsim 

-^0 

2.509 

1.901 

5.997 

5.694 

12.499 11.749 

EMCt 

5.002 

3.993 

11.565 

10.775 

21.748 21.082 

j-^approx 

^0 

6.224 

5.705 

16.535 

15.831 

37.159 36.564 


schematic MCT, which is known to overestimate the coupling constant A. However this 
analysis not only explains the behaviour of MCT at high temperature, it also throws some 
light in the origin of its breakdown at low temperature. Usually the breakdown of MCT at 
low temperature has been attributed to the neglect of higher order correlation functions^i^. 
This present analysis predicts that the stronger temperature dependence of the vertex might 
be partially responsible for the breakdown of MCT even at low temperature. 

B. The Adam Gibbs Relation and MCT 

We have shown that the relaxation time ,r, over a temperature regime (10“^ < (^ — 1) < 
10*^) follows both the AG relation and MCT power law behaviour. We also find the avoided 
divergence obeserved in the configurational entropy plot (Figl2]) arises from the vanishing of 
the pair configurational entropy {Sc 2 )- For all the systems studied here, we find Tk 2 — Tc 
(Table [I]), thus we can rewrite EqJTHl as. 

- 1 ) =; KT2{f, - 1 ) ( 31 ) 

XK2 J-c 

We note that although Tk 2 — Tc and the MCT framework which predicts the power law 
behaviour is developed at the two body level, the AG relation with Sc 2 alone cannot predict 
the MCT power law behaviour, and the RMPE ,A5', plays an important role in predicting 
it. We also show that indeed there is a relation between MCT critical exponent 7 , Adam 
Gibbs coefficient A, the pair thermodynamic fragility /Ft 2 - 

As shown earlier in EqJT^ conhgurational entropy can be written in terms of pair con¬ 
figurational entropy and RMPE. Thus we can write. 
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A 


A 


A 


1 


(32) 


TS^ TS^2 + TAS Kr2 1(^-1) + ^] 


where we have used EqUS] and Eq|3l] to write the hrst and second equality respectively. 

We hnd that although TAS is system dependent (FiglHK), except for WCA system at 
p = 1.2 the function shows a master plot when plotted against — 1) (Fig| 8 )D). Note 
that although the value of is small, it is not negligible. 

The master plot of can be htted to a straight line, = 0.26 — 0.35(^ — 1). Next 
we show that a function when plotted against ln(^ — 1 ) shows linearity in the 

whole regime of (10“^ < — 1) < 10*^) only when /(T) is non-negligible positive quantity 

(FiglHl::). Note that in FiglHl:: when /(T) = 0 (which implies AS* = 0 in Eq l3^ the function 
diverges strongly. This shows that the AG relation at two body level cannot predict the 
MCT power law behaviour. 

The analysis further shows that to obtain a correct estimation of the MCT power law 
exponent 7 (slope of the plot), f{T) needs to obey the following temperature dependence, 
/(T) = = 0.26 — 0.35(^ — 1). The two functions and (^ — 1) show opposite 

trends, the former increases whereas the later decreases with temperature. Therefore a 
crossover between these two functions is observed in this regime and around MCT transition 
temperature, S> (^ — 1 ) and conhgurational entropy and the relaxation time are 

determined primarily by many body contributions. 

From FigJHli we hnd in the temperature regime (10“^ < (^ — 1) < 10°) Eql32] can be 
re-written as. 


A A 1 mA T 


(33) 


where ‘m’ is the slope obtained from FigIHl(d) and given in Table IIIIl Since r is found to 
follow AG relation we can write. 


T ~ “P(2^) ~ (jT - 1)*" 


(34) 


Comparing Eq|l4] and EqlM] we can write. 


mA 

Kt2 


rsj 


7 


(35) 
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FIG. 8 . (a) TAS is plotted as a function of — 1) and it shows a strong system dependence, 

(b) vs — 1) showing a master plot for all the systems except for WCA system at p = 1.2. 
Dotted line is guide to eye. (c) ((^ — 1) + /(T))“^ plotted against — 1) by varying f{T). Only 
for non negligible positive values of f{T), linearity is found in the regime 0.1 to 1.0 of — 1). 
To obtain a correct estimation of the MCT power law exponent 7 (slope of the plot), f(T) needs to 
be temperature dependent (green dashed line), (d) [jr — ^ ta$ ] vs — 1 ) shows a master plot 
for all the systems except for WCA system at p = 1.2. ‘m’ is the slope of the linear plot which is 
tabulated in Table \IB 


where m ~ 1 for all the systems except for the WCA system at p = 1.2. Thus we show 
that the MCT scaling parameter, 7 is related to the AG parameter, A and the pair ther¬ 
modynamic fragility of S'c 2 , Kt 2 - We have tabulated the 7 values in Table HVl which shows 
the above relation holds. The deviation of slope value (‘m’) from unity for WCA system at 
p = 1.2 may have some connection to its breakdown of density-temperature scaling which 
needs to be investigated in future. 
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TABLE III. The slope of the linear plot o/[--y— ^ jAg ] vs (^ — 1) in the region (10 


-1 


10°) (FigWl)- 


p 

m{LJ) 

m{WCA) 

1.2 

0.987 

0.695 

1.4 

1.029 

0.888 

1.6 

1.004 

1.000 



1 )< 


TABLE IV. mA/KT 2 o,nd 7 for LJ and WCA system. As predicted by Eaf3^ mA!K to value is 
similar to 7 value obtained from free fitting (Fig[l\ ) for most of the systems. 



p = 1.2 

p= 1.4 

p = 1.6 


mA/KT2 

7 

mA/KT2 

7 

mA/KT2 7 

LJ 

2.322 

2.229 

2.474 

2.385 

2.352 2.299 

WCA 

2.932 

2.243 

2.852 

2.289 

2.579 2.304 


The MCT critical exponent (7) is known to be density-temperature independent^. Inter¬ 
estingly we also find that although both AG coefficient (A) and pair thermodynamic fragility 
{Kt 2 ) are strongly dependent on density and temperature (Table IV)) . but their ratio, which 
is related to 7 (Eq l3^ . is density-temperature independent (Table [TVT) . 


VII. CONCLUSION 

In this work we show that in a certain region (10“^ < — 1) < 10°) the relaxation 

time follows both the AG relation and MGT power law behaviour. We also find that the 
MGT divergence temperatures coincide with the temperature where pair configurational 
entropy goes to zero for all the systems studied here. AG relation is based on activated 
dynamics, whereas MGT is mean field theory which at the two body level does not address 
any activated dynamics. Also the microscopic MGT does not have any apparent connection 
to entropy. Thus to understand the above mentioned observations we explore the connection 
between mode coupling theory and entropy and discuss different predictions of MGT in the 
light of entropy. 

In this article we show that the MGT vertex for the structural relaxation time under 
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TABLE V. The Adam Gibbs coefficient ‘A’, as obtained from the linear fit of t vs 1/TSc plot, and 
pair thermodynamic fragility Kt 2 , cls obtained from the linear fit ofTSc 2 vs T/Tk 2 plot for both the 
systems at different densities is tabulated below. The data shows that both are strongly dependent 
on density. 


p 

A{LJ) 

A{WCA) 

Kt2 {LJ) 

Kt2 {WCA) 

1.2 

1.87 

1.89 

0.795 

0.483 

1.4 

3.57 

4.37 

1.555 

1.358 

1.6 

6.96 

7.57 

2.971 

2.936 


certain approximations can be related to the pair excess entropy. Higher order MCT calcn- 
lations in the schematic MCT framework can relate the relaxation time to the exponential 
of this vertex. Thus the MCT can provide a microscopic derivation of the phenomenolog¬ 
ical Rosenfeld theory. Our analysis shows that the Rosenfeld parameter is related to the 
static structure factor S{k). The temperature dependence of S{k) leads to the temperature 
dependence of Rosenfeld parameter ‘K’, thus explaining the earlier observation of the non¬ 
uniqueness of the Rosenfeld exponent^'^^. The analysis of the vertex reveals that quantity 
which contributes to the vertex, S 2 apporx has a much lower value and stronger temperature 
dependence as compared to the excess entropy , Sex- If we assume the Rosenfeld scaling to 
be valid and replace Sex by S 2 approx, the predicted relaxation time shows similar character¬ 
istics as the MCT relaxation time. Thus the study reveals that the larger value of tmct 
and its higher activation energy is related to the value and temperature dependence of the 
vertex. This analysis further reveals that the breakdown of MCT at low temperature might 
be partially related to the strong temperature dependence of the vertex. 

As mentioned earlier the AG theory which is based on activation dynamics can com¬ 
pletely describe the mode coupling theory (MCT) power law behavior in the region where 
the latter is found to be valid. Since the configurational entropy has a finite value at the 
MCT transition temperature, Tc, the AG relation is not expected to predict any avoided 
transition in this regime. Our study reveals that although Sc is finite, Sc 2 vanishes at Tk 2 
(where Tk 2 = Tc), thus being responsible for the divergence like behavior. However we 
show that the pair configurational entropy although predicts the correct MCT transition 
temperature it by itself cannot predict the MCT power law behaviour. The residual mul- 
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tiparticle entropy (RMPE) plays an important role in providing the correct temperature 
dependence of relaxation time. We also obtain a connection between the AG coefficient (A), 
pair thermodynamic fragility {Kt 2 and MCT critical exponent (7) and found although first 
two quantities are dependent on density and temperature, their ratio which is related to 7, 
is density-temperature independent . 

Note that although the absolute value of AS is in the similar range both at high and low 
temperature regimes, in the high temperature regime it plays a minor role in determining the 
dynamics, whereas its role at low temperature becomes central as we approach the avoided 
transition. This small positive value of AS playing an important role in predicting the MCT 
power law behaviour is similar to the prediction of unihed theory^. In the unihed theory 
it was shown that in a certain temperature regime many body activated dynamics plays a 
hidden but central role in predicting the MCT like behaviour of the total relaxation time. 
Although apparently the MCT does not depend on the properties of landscape, the saddles 
in the landscape have been found to disappear at . Here we show that Sc 2 also vanishes 

at Tc- Thus there may be a connection between pair configurational entropy and saddles. It 
will be also interesting to understand the independent role of pair conhgurational entropy 
and RMPE in the landscape picture. These are important open questions to be addressed 
in the future work. 
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